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2^ ■ Abstract 

We introduce a general method for transforming the equations of motion following 
Qh. from a Das-Jevicki-Sakita Hamiltonian, with boundary conditions, into a boundary 

,^ ■ value problem in one-dimensional quantum mechanics. For the particular case of a one- 

dimensional chain of interacting N x N Hermitean matrices, the corresponding large 
^ ' A^ boundary value problem is mapped into a linear Fredholm equation with Hilbert- 

5r , Schmidt type kernel. The equivalence of this kernel, in special cases, to a second order 

differential operator allows us recover all previously known explicit solutions for the 
matrix eigenvalues. In the general case, the distribution of eigenvalues is formally derived 
through a series of saddle-point approximations. The critical behaviour of the system, 
including a previously observed Kosterlitz-Thouless transition, is interpreted in terms of 
the stationary points. In particular we show that a previously conjectured infinite series 
of sub-leading critical points are due to expansion about unstable stationary points and 
consequently not realized. 



1 Introduction 

The utility of studying the statistical mechanics of systems which can be encoded in terms 
of matrix variables has been evident over the past twenty years. From the generalization of 
Wigner and Dyson's work by Brezin, Itzykson, Parisi and Zuber [^ to the use of matrix models 
in low dimensional gravity (see [^ for reviews), to modern implementations of matrix strings 
0, matrix models have proven to be a convenient framework for calculation. Remarkably, 
throughout this period the dynamics of matrix models have been considered only in the 'sta- 
tionary' case where solutions of the equations of motion are assumed to be time-independent. 
There arise circumstances though where the time evolution of a system governed by an action 
over matrix variables is relevant, especially if there are definite initial and final conditions on 
the state of the system. Examples include two-dimensional Yang-Mills theory on the cylinder 
and sphere 0, |^, matrix model approaches to the continuous Hirota equation [Q and lattice 
QCD 1^ and, the case we will focus on here, linear chains of interacting Hermitean matrices. 
It is our goal to give a coherent framework in which in to study this class of matrix model 
problems. 

When an action over matrix variables can effectively be reduced to the eigenvalues, the 
treatment of the eigenvalues as a single, collective coordinate [^ provides a convenient view- 
point for determining the Hamiltonian of the system. For the matrix actions we will consider 
here the particular form is known as the Das-Jevicki-Sakita Hamiltonian 0. The associated 
equations of motion are the non-linear Euler equations of fluid dynamics in one dimension 
and with the addition of appropriate boundary conditions present a difficult mathematical 
problem. Here we will show how these boundary value problems can be interpreted in the 
familiar framework of the semi-classical limit of one- dimensional quantum mechanics. While 
in general solutions are still difficult to obtain, we will focus our attention on the case of a 
one-dimensional chain of interacting N x N Hermitean matrices which was shown to be a 
Das-Jevicki-Sakita-type system by Matytsin |T^ in the limit of large A^. In this situation a 
restatement of the equations of motion in terms of quantum mechanics allows one to consider 
a linear problem which is most conveniently stated as a Fredholm integral equation. Formal 
solution of this equation in terms of iterations of the integral kernel, which can be evaluated 
by saddle-point methods, will provide a framework for detailing the nature and position of 
critical points which arise in the system. 

An interesting system from a statistical mechanical point of view, the one-dimensional 
matrix chain was discussed previously as a model for string theory with a discrete target 
space IT^. In addition to this direct application, the matrix chain is related by duality 



transformation to string theory on a circle with radius R = 1/a |T^, |T^, string theory at 
a finite temperature [ffH] and to 0(2) sigma models coupled to two dimensional quantum 



gravity [|T6i p!7| . The model itself consists of Hermitean matrices (px on the sites {x} of a 



one-dimensional lattice interacting with nearest neighbours. The action of the system is 

X X ^'^ 

where a is the lattice spacing and V is the potential at each site. Assuming the site potentials 
V{(j)x) to be analytic in (px, the action can be reduced to the eigenvalues of the matrix variables 
by SU{N) gauge rotations A and the known result for the correlator [|ll|] 

I[<Px, <Px+l] =JdA e-^TrA0.At0.+, ^^2) 

Consequently, in the limit of large A^, the calculation of the partition function reduces to 
determining solutions of equations of motion for the eigenvalues. Despite its relative simplicity 
and importance in low-dimensional string theory and quantum gravity, explicit solutions of 
the matrix chain are only known for a few particular choices of potential V. In Section 
Four we will show these cases arise when the equations of motion are equivalent to certain 
Sturm-Liouville eigenvalue problems. 

This model exhibits a remarkable range of behaviour as one adjusts the lattice spacing, a. 
For instance, in the limit of vanishing a one recovers the canonical example of a Das-Jevicki- 
Sakita system, matrix quantum mechanics (restricted to a singlet sector) 

Smqm = jdx Tr[i02(x) + V{ct>{x))] (1.3) 

In fact, the methods we will introduce in the next Section and use to study the matrix chain 
find their most natural form when applied to such matrix mechanical systems on a finite 
interval x G [0,T]. 

More interestingly, fixing the potential V to satisfy V'ifi) = —4, at critical lattice spacing 



a = 1, the matrix chain was argued to undergo a Kosterlitz-Thouless phase transition 113 
The physics of this phase transition are most transparent in the dual picture of string theory 
on the circle. Here strings wrapping around the non-contractible target space play the role 
of vortices which are liberated when the radius R ~ 1/a becomes sufficiently small. In terms 
of a discrete lattice, at sufficiently small lattice spacing a < 1 the universal behaviour of 
the system is that of the continuous model (|1.3| ) due to communication between lattice sites. 
When the spacing is larger than unity this communication is broken and the system can be 
expected to factorize into a product of single D = matrix models, which in fact it does as 
a — s> oo. The main objective of studying the matrix chain is to understand precisely how this 



change of character from D = 1 matrix quantum mechanics to isolated D = matrix models 
takes place at a = 1. 



In addition to a = and a = 1, it was argued in |T6| and [^ that at a discrete infinity of 
lattice spacings non-analytic behaviour arises in the distribution of eigenvalues of the matrix 
variables (px- Considering iterative approximation of a functional solution to the equations of 
motion, it was argued these critical points arise when a = sin {piT/2q) where p < q are each 
positive integers. In Section Five, where we develop the general solution for the distribution 
of eigenvalues as a series of saddle-point integrations, we will show that this behaviour is not 
realized in the matrix chain. In fact it will be shown that such critical behaviour arises when 
unstable stationary points are mistakenly considered relevant in saddle-point approximations. 
Moreover, in constructing the general solution it will become clear that a countable infinity of 
unstable stationary points appear in intermediate steps each of which can lead to misleading 
results. 

In addition to the details just outlined, we recall the Kazakov-Migdal (lattice) model 
of induced QCD as an interesting application of the techniques developed here. Since the 
D-dimensional model is effectively reducible to the one-dimensional case, the non-trivial be- 
haviour of the matrix chain demonstrated here may help to overcome some of the apparent 
deficiencies of that model of gauge interaction. 

2 Collective field theory to quantum mechanics 

To begin we will demonstrate a non-linear change of variable which allows the equations of 
motion derived from a Das-Jevicki-Sakita Hamiltonian [H 



n[p,n] = ydxpix)l(^^) -^p\x)\ + ldxp{x)vix,t) (2.1 



to be determined as a temporal boundary problem in one- dimensional quantum mechanics. 
In ( |2.1| ), p(A) is a probability distribution function which returns the fraction of eigenvalues 
of the matrix variables (p^ with value A. In the present context of collective field theory, p 
plays the role of a coordinate and 11 is the canonically conjugate momentum. Here we have 
included a potential V which is most generally time- dependent. Independent of time, ( |2.1|) 
describes the evolution of the matrix quantum mechanics ( |1.3| ) with V = V. 

The equations of motion following from this Hamiltonian are given by the functional 



variations with respect to the canonical variables 



d 571 d d (6n\ 

where, for convenience, we have defined the 'velocity' field 

v{x,t) = — (2.3) 

Explicitly we see that the equations of motion are of the form of the Euler equations for a 
one-dimensional inviscid fiuid experiencing a time-dependent force derived from V 

| + ^(P")=0 (2.4) 

dv 1 fdv^ 2<^P^\ "^ 

dt 2 \dx dx J dx 

As in all evolution problems we need to specify boundary conditions to supplement the 
equations of motion. In the case of matrix quantum mechanics the classical treatment of |l| 
assumed free boundary conditions for the evolution of the system over infinite time. In this 
case one is led to a study of the ground state of the system and the equations of motion ( |2.4| ) 
are solved for the stationary case of f = 0. In general though we are interested in arbitrary 
boundary conditions for the evolution of the system through some time T. To be definite we 
will specify the generic conditions on the fields v and p at time t = and t = T 

p{x,0) = poix) , v{x,0) =vo{x) (2.5) 

p{x,T) = Prix) , v{x,T)=vt{x) 

We would now like to rewrite the equations of motion in terms of a single quasi-linear first 
order differential equation. Introducing the complex function 

f{x,t)=v{x,t)+i7!'p{x,t) (2.6) 



it is easy to see that (|2.4|) leads to 



dt dx dx 

For vanishing potential and real / this is the Hopf equation which is the prototypical one di- 
mensional model of an equation which admits wave solutions complete with shocks and other 



inherently non-linear behaviour |T^. We would like to solve this equation with the boundary 
conditions given in ( p.5|) . While the solution of the (complex) Hopf equation is known explic- 
itly |T^ in terms of initial data, typically we are interested in cases where boundary data is 



supplied. In these situations the exphcit solution leads to an ill-posed, inverse-type problem 
for the boundary conditions. Consequently we would like to find a more natural presentation 
of the problem. To this end we will consider a generalized version of (p.7|) and introduce a 



dispersive term [^ which will serve to smooth out the solutions and allow one to calculate 
the details of the solution / in terms of boundary conditions. Adding such a term leads to a 
complex version of the forced Burgers' equation 



ih 






^ + /^ + Av(x t) 
dt dx dx 



{21 



Obviously, in the limit as /i -^ the forced Hopf equation (|2.7| ) is recovered. The useful- 
ness of this addition to the problem is that there exists a well-known (to mathematicians) 
change of variable which will reduce ( p.8| ) to a second order, linear differential equation. This 
transformation is known as the Cole-Hopf transform and is given 



d 

f(x,t) = -2ih—log'4j{x,t) 



(2.9) 



As is readily verified, this transformation leads (up to an irrelevant additive constant of 
integration) to the Schrodinger equation 



Hij 



-h' 



(92 V(x,t) 



+ 



ij 






(2.10) 



dx^ ' 2 

Consequently we see that the time evolution of the collective field Hamiltonian ( |2.1|) is given 
by the semi-classical {h — > 0) limit of one-dimensional quantum mechanics with the same 
potential, up to a factor. 

Restricting to time independent potential^ V{x,t) = V(x), the boundary conditions (|2.5|) 
enter into the picture as the start and end points of evolution by H through time T 



^(x,T)= e-^^^/V(a;, 0) 
where the wavefunction is defined in terms of the eigenvalue density and velocity 



ip{x,t) = exp 



-^fdvPiv,t) + ^fdvviv,t)^ 



(2.11) 



(2.12) 



Given the Green's function G{x,u\T) for the Schrodinger problem it may be convenient 
to cast ( p.ll| ) as an integral equation 



exp 



-^ f dv pAv) + ^ f dr^ vriv) 

/ du G{x,u\T) exp 



(2.13) 
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2h 



dv Poiv) + 



2h 



dr] voiv) 



^This restriction is not essential. In the general case the time evolution of the system is given by a Dyson 
series expansion. 



In applications the implementation of boundary conditions and Green's function varies de- 
pending on the situation. For example, in investigations of the continuous Hirota equation ||^ 
one is interested in the case of time-dependent potentials V{x, t) leading to non-trivial Green's 
functions. Here the initial and final velocities are specified and the initial and final eigenvalue 
distributions are identified (po = Pt) so that ( |2.13| ) takes on the form of a non-linear integral 



equation for po- 

Another, more familiar application is that of two-dimensional Yang-Mills theory 0, ^ |T9|. 
In this case the coordinates (x, u) in ( |2.13| ) are replaced by their periodic counter-parts since 
the matrix variables are unitary rather than Hermitean. The Green's function is that of the 
heat equation on a circle and the boundary conditions on p and v specify the topology of the 
two-dimensional space-time. In general these cases again lead to non-linear integral equations 
for the unknown fields but as we will now see, there are situations where the general integral 
equation ( |2.13| ) leads to a tractable, linear equation. 

3 An application: the matrix chain 

An application of the preceding formalism, which we will be investigating for the remainder of 
this Paper is that of the one- dimensional matrix chain. Instead of ( |1.1|) , it will be convenient 
to consider a different form of the action 

Su = Yl Trf/(</).) - Y. Trc^x'/'x+i (3-1) 

X X 

where we have rescaled the matrix field (px and the new site potential U is related to V by 

t/(0j = a\/(v^0,)+02 (3.2) 

The evolution of eigenvalues in this system, in the limit of large A^, was shown by Matytsin 



T0| to be given by the Das-Jevicki-Sakita Hamiltonian ( |2.1| ) with vanishing potential V. Here 



the boundary conditions (|2.5|), which play a crucial role, are given by 



Po(x) = p(x,0) = p(a;, 1) (3.3) 

v{x,0) = —v{x,l) = -U'{x)—x 

The objective is to solve for the initial and final eigenvalue density po in terms of the potential 
U. 



With V = 0, and e = ifi, the Schrodinger equation (|2.1CI|) reduces to the linear heat 



equation. Substituting the well known Green's function in ( |2.13| ) and implementing the 

6 



boundary conditions ( p.3| ) , we find a linear integral equation of second Fredholm type for the 
initial eigenvalue density po{x) = p{x, 0) in terms of the effective potential U{x) 

^(^) = -7^ due ^^ ^ 4^ x{u) (3.4) 

v47re "' 



where 

in 



X{x) = exp 



2 , dr] po{7]) 



(3.5) 



and we have introduced the constant A which does not affect the asymptotic solution but 
is convenient for later discussion. In this form ( |3.4| ) has the natural interpretation of an 
eigen-problem with eigenvector x and associated eigenvalue A. It is interesting to note that 



the same integral equation arose in the context of matrix chains previously [0, |T3| . In these 
instances the authors were seeking the lowest lying eigenvalues A in order to calculate the 
partition function. Here we are most concerned with the asymptotic (e — > 0) form of the 
eigenfunctions x ^^ order to recover po- 

Hence we have re-stated the differential form of the boundary value problem (eqns. (|2.4|) 
and (|2.5|) ) for the matrix chain as a well-posed, linear problem for the eigenvalue density po- 
Moreover, the integral form ( |3.4| ) allows one to write equations for the eigenvalue density 
immediately. Formally, in the limit of vanishing e, the integral in ( ^.41 ) can be calculated 
by saddle-point methods. This calculation is particularly straightforward when only a single 
stationary point, Us contributes to the integral ( p.4|) . This stationary point is determined by 
the extrema condition 

X = U'{us)/2 + mpoius) (3.6) 

Likewise, at this stationary point po can be formally evaluated using its relation to x from 
equation ( |3.5| ) 

Us = U\x)/2-mpQ{x) (3.7) 

Labeling x by 6*+ in equation ( p.6| ) and Ug by G_ in equation (|3.7|), we see the saddle-point 
evaluation of the integral equation p.4| ) leads to Matytsin's solution [|TU] of the boundary 
value problem 

x = G+{G^{x)) = G^{G+{x)) (3.8) 

This form of the solution of the problem is well-suited to detailed local analysis as demon- 



strated in ll^, |T^. However, for a global analysis which can accommodate the existence of 
multiple stationary points, we will take the (well-studied) integral form ( p.4| ) over the func- 
tional equation ( |3.8| ). Before considering the particular characteristics of solutions though, 
we will utilize the interpretation of integral operators as self-adjoint operators to recover all 
previously known explicit solutions po of the matrix chain. 

7 



4 Integral and differential operators 

An analysis of an integral equation centers on the properties of the kernel K. In our case, the 
kernel of (|3.4|) is 

r^/ \ 1 ^^ Uix) + U(u) 

K{x,u) = ea^ 4f (4.1) 

V47re 

This positive definite kernel is real, symmetric and hence Hermitean. In the following we will 

assume the kernel to have many convenient properties which follow from the boundedness of 

all n-fold iterations 

/ dxi ■ ■ ■ dXn~lK{xi,X2) ■ ■ ■ K{Xn-l, x„) < oo (4.2) 

Consequently we will only consider potentials U{x) which increase at infinity faster than x^. 
Under such conditions K defines a compact operator and represents a self-adjoint operator 
on the space of square-integrable functions. From classical operator theory [Q there is a 
countable infinity of orthogonal eigenf unctions {Xn{x)} which form a complete set on this 
space. It is our goal to obtain, as a function of the site potential U, the asymptotic form of 
these eigenf unctions which determine the large A^ solution of the eigenvalue distribution po 
through (p.5|). 

There are a number of cases in which the kernel K is equivalent to familiar self-adjoint 
operators in physics. In fact only in these instances are explicit solutions of the large N 
matrix chain known. In each case the eigenvalue distribution po is the solution of a quadratic 
equation. Here we will show that such situations arise by considering the e — > limit of 
second order self-adjoint differential operators C for which C and K share a complete set of 
eigenfunctions and define the same action on the Hilbert space. 

In order to determine the conditions that C must satisfy to be equivalent to the integral 
kernel we act on Xn as defined by the integral equation 



EnXn{x) = C^Xn{x) = J du C^K{x,u)Xn{u) (4.3) 

Likewise, acting under the integral sign and using the self-adjoint property of £ 

EnXn{x) = I du K{x, u)CuXn{u) = j du Xn{u)^uK{x, u) (4.4) 

Subtracting ([4.4|) from ([4.3|) we find that C will be equivalent to the integral kernel K, if and 



only if we satisfy the natural commutativity condition 

{C,-Cu)K{x,u) = (4.5) 



It follows from the symmetry of the kernel that the action of L generates symmetric functions 
of X and u which will cancel out in 



CxK{x, u) = h{x + u, xu)K{x, u) ='^ amn{x + u)"^{xu)"^K{x, u) (4.6) 

mn 

where h is an arbitrary function of two variables with Taylor expansion coefficients {amn\- 

This construction is valid for any differential operator but for simplicity we will consider 
only second order differential operators of Sturm-Liouville form 

^^ = ~^''^ X'^^^'^iij + '''^''^ ^^-^^ 

Now it is a simple matter to find the potentials U and associated coefficients 02 and Oq which 
are consistent with the symmetry condition (|4.6| ). Since a second order C can only produce 
quadratic polynomials in u, (|4.6|) leads to the equation 



-a2{x){u-U\x)/2f-2e{a2{x){u-U\x)/2))' + a^{x)= ^ a^n{x + u)"^ {xuY (4.8) 

m+n<2 

Expanding each side of the second equality in powers of u and equating coefficients gives, to 
leading order in e, the result 

^ji, ^ _ aiix^ + (2a2o + Qoi)^:^ + Qio / . qx 

CiQ2^ + OllX + (320 

02(0;) = — (ao2a;^ + aiix + 020) 

,2 {aux^ + (aoi + 2a2o)x + aiof 

ao[x) = aoo + aiox + a2oX — ■ r 

4(ao2a:"^ + aiix + a2o) 

Having determined the coefficients 02 and oq, we can solve for the explicit form of the asso- 
ciated eigenvalue distribution. The differential form of the integral equation (|3.4|) is 



C.xix) = Ex{x) (4.10) 

where the eigenvalue E is related to that of the integral form by A = e^/^. In the limit of van- 
ishing e, the solution of this differential equation is given simply by the WKB approximation 
which implies the replacement of derivatives with 'momenta' 

With this substitution, ( [4.10| ) reduces to a quadratic equation in the eigenvalue distribution 
which is solved by 



Poia;J 



E - aoix) ^^^^2) 



n\ a2 



X) 



4:{ao2x'^ + aux + a2o)(aoo + ctioa:; + ^202;^ - E) - (anx^ + [am + 2a2o)a; + aio)^ 



2TT{ao2x'^ + aux + 020) 



Since po is to be interpreted as a probability distribution describing the eigenvalues of the 
matrix variables, it must be properly normalized. This is carried out by integrating po over 
the positive support of (ag — E)jai and fixing the undetermined constant E to satisfy 



dxpo{x) 



(4.13) 



Integrating (|4.9| ), we find the most general potential of the matrix chain which is consistent 
with a second order, self-adjoint differential operator 



U{x) 



flu 2aQ2aio — aoiflo2«ii + o-u — 2002^11020 
X — . arctan 



Q02 



Q01Q02 ~ Qii 

2^02 



O-n 



log 



'02 V ^'^02^^20 
«20 + OllX + aQ2X 



On + 2ao2a^ 
4ao2'^2o ~ c^ii 



(4.14) 



Contained in this general form are a number of previously examined examples of the matrix 
chain. First, setting an = and aio = recovers the symmetric double Penner-type potential 



first considered by Matytsin [|T0[. Also the asymmetric double Penner considered in connection 
with the Kazakov-Migdal model of induced QCD PT| can be recovered (see 



for recent 



work). In each of these cases the eigenfunctions x ^i-re related to prolate spheroidal wave 
functions. Additionally, with 020 = 1, the limit ao2 -^ recovers the quadratic potential case 
solved by Gross p3[ and Makeenko [Q. With U' = (aoi + 2)x = 2rn?x, C is the Schrodinger 
operator for the simple harmonic oscillator and the general expression for the distribution of 



eigenvalues ( [4.12 ), properly normalized, reduces to 



TT 



Po{x) = -V2Vm4 - 1 



m^ 



l)x2 



(4.15) 



In each of these cases the analytic structure of the eigenvalue distribution (|4.12|) is what one 
would expect of a D = one-matrix model with a particular potential. This correspondence 
was previously suggested for arbitrary potential in the matrix chain, but such a simple solution 
of the integral equation ( |3.4| ) is not possible. In fact with further computation it can be 
argued that higher, finite, order self-adjoint differential operators are equivalent to the integral 
kernel K only if they are functions of the second order operator C^ we have constructed 
above. Consequently, the large N eigenvalue distribution po of the matrix chain is the solution 



of a polynomial equation only for potentials of the form ( |4.14| ). This suggests that for a 
generic potential the solution of the matrix chain problem follows from the semi-classical 
approximation to a pseudo-differential operator as suggested by naively extending the results 



for finite chains of matrices ||2^. We will not pursue this line of reasoning here but return 



in the next Section to the integral form (O) and give a method for obtaining the general 
solution and describe its unique features. 
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5 General solution and classification of critical points 

A phase transition in a system is characterized by a change in the analytic structure of 
observables as functions of external parameters. In standard matrix models this is commonly 
taken to include a change in analytic structure of the distribution of eigenvalue^. In the 
present case we have shown that the distribution of eigenvalues for a linear chain of Hermitean 
matrices is determined completely by a linear integral equation ( |3.4|) involving the kernel K 
( |4.1| ). In this Section we will go further and solve for the asymptotic eigenf unctions x as a 
series of iterations of K. Given that iterations of K are calculable by saddle-point methods 
and the relation (^.51), we can effectively solve for the distribution of eigenvalues po- In 
particular we will be able to determine the critical structure of the matrix chain by analyzing 
changes in the analytic structure of the kernel K and its iterations. 

We begin by slightly generalizing the problem of finding the eigenfunctions x i^i order to 
make use of well-known techniques [2^ . Instead of the integral equation (|3.4|) let us consider 



for a moment an inhomogeneous version with c some real constant 

x(x) = c + A j du K{x, u)xiu) = c + XKx (5.1) 



Since there are no constant eigenfunctions of the kernel K, the unique solution of ( |5.1| ) is 
given by the Neumann series which builds the solution up through an iterative procedure 

Xo = c (5.2) 

Xi = c + XKxo 

Xn = C+ XKXn-1 

Taking the limit of this process generates the solution 

X{x) =il + Y.r I du K"(x, u)\c = R{x; X)c (5.3) 

where R is commonly referred to as the resolvent kernel and i^" is given by the convolution 
of n kernels 

fC"(x,M)= dzi-- ■dzn-iK{x,zi)K{zi,Z2) ■■ ■K{zn^i,u) (5.4) 



^This assumption of a physical phase transition following from a change in the analytic structure of 
the eigenvalue distribution does not always hold |26|. Nevertheless, in the absence of calculating physical 
observables, we will take a change in structure of the eigenvalue distribution to be a strong hint of a physical 
transition. 
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The convergence of the resolvent kernel is guaranteed by the boundedness of -ft"" and the 
Hilbert-Schmidt theory of integral equations |2^. In order to recover the homogeneous solu- 
tion we should set c = which leads to a trivial solution for x- Fortunately for us only the c 
independent logarithmic derivative of x is required to recover the eigenvalue distribution po 
of the matrix chain. From the definition of Xi (|3-5|) and ([5 .31) we have the result 



Pq{x) 



li"?i — 1~ l°g ^(^) = li^ — T 1°S R{x] A) 

«^0 TT dX <:^0 TT dX 



(5.5) 



Hence we can extract the eigenvalue distribution from the resolvent kernel R, which is 
expressible solely in terms of iterations of the kernel K as defined in (|5.3| ). This is an important 
observation since the iterated kernels are calculable in the limit of vanishing e. Explicitly, the 
^th iiQYdXed kernel can be written for n > 2 



, , 1 U(x) + U{u) f _„, 



K''(x,u) 



(5.6) 



where the 'action' S is given by 

n- 



n-1 



n-2 



(5.7) 



1=1 



i=l 



Evaluating the integral over z = {zi, . . . , Zn-i} is a straightforward application of saddle- 
point methods in the limit of vanishing e. In this limit the dominant contributions to the 
integral come from the vicinity of the stationary points of the function S 



dS_ 
dzi 







(5.^ 



In general there are many such stationary points, which we will label by a. Expanding S up 
to quadratic order about the stationary points, the resulting Gaussian integrations give the 
approximation to the iterated kernel 



K"+Vx,m) 



V4: 



1 UM + U(u) (-l^niS^) 

e 4e 2^ , = e 



S{z")/2e 



Tce 



detHniz" 



(5.9) 



where A„(z") is the number of negative eigenvalues of the nx n Hessian Tin of second deriva- 
tives of S evaluated at the stationary point 



njz' 



d^S 



dzidzj 



I U"{z^) -1 ■ ■ ■ \ 
-1 U"{z^) -1 

V ••• 0-1 U"{zZ) ) 



(5.10) 
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While we have counted contributions from all stationary points of S in ( p.9| ), the only relevant 
ones are those with minimum 'action'. The contributions of other, irrelevant, stationary points 
are exponentially suppressed in the limit of vanishing e. 

Unfortunately it is difficult to calculate the iterated kernels and the resolvent in closed 
form to obtain explicit results for the eigenvalue distribution. It is possible though to extract 
enough qualitative information from ( |5.9|) to outline a one-dimensional phase diagram of the 
matrix chain. To be definite we will consider a quartic potential 

U{(j)) = (j)^ + m^(f)'^ (5.11) 

This simple form is convenient since it allows a clear view of the intrinsic behaviour of the 
matrix chain which depends on the strength of the quadratic term. 

Starting with m^ large and positive we expect that the eigenvalue distribution will be 
localized near the origin where the quadratic term of the potential is dominant. Consequently, 
the solution (|4.15|) for the pure quadratic potential is a good approximation to the true 



solution and is smooth deformation of it. Making contact with the saddle-point evaluation 
of the iterated kernels ( ^.9|) it is easy to convince oneself that for large m^ there is only one 
stationary point of the action ( ^.7|) and it is located near the origin z = 0. As m'^ is decreased, 
the solution ( ^.15| ) of the pure quadratic potential becomes unstable at m'^ = 1 and in fact 



this instability is present in the full solution. Information about the stability of stationary 
points is contained in the eigenvalues fir of the Hessian 7is{z), which evaluated at the origin 
in z-space give 

/ir = 2m^ — 2cos , r = l,...,s (5.12) 

Consequently we see that when evaluating arbitrarily high iterations of the kernel, K^, the 
Hessian at the origin will being to develop negative eigenvalues at rri^ = 1 and this stationary 
point becomes unstable. Since S is bounded from below for our choice of potential, there are 
other stationary points where S has a lower magnitude and it is these which will dominate 
the integral (|0| ) for K'^. 

This degeneration at m^ = 1 signals a phase transition in the matrix chain and the physical 
meaning of it can be found by returning to the original form of the matrix chain action ( |1 . 1|) . 
The potential there, V is related to the current U by the one dimensional lattice spacing a 
through 

U{(l)) = aV{y/E(f)) + (f)'^ (5.13) 

Taking two derivatives with respect to and fixing (without loss of generality) V"{0) = —4 



to agree with the conventions of |1T6|], we find 

m'^ = l-2a^ (5.14) 
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Hence, the instability at m'^ = 1 corresponds to the hmit of vanishing lattice spacing and we 
are observing a phase transition to D = 1 matrix quantum mechanics (i.e. ( |1.3| )) from the 
m^ > 1 phase where the solutions are qualitatively what one would expect from a. D = 
single matrix model. In fact this observation can be made more concrete by returning to the 
integral equation (p.4|) but in terms of the potential V and lattice spacing a, 



e-ir^. Mv) = ^fdu e^'^-^'^'^'i^ e'^f'^ ^^^^^ (5.15) 

It is a straightforward calculation to show that for vanishing e, and to leading order in van- 
ishing a, ( ^.15| ) is solved by 

poix) = -Je-V{x) (5.16) 

which is the well-known solution for the ground state of matrix quantum mechanics |^ . 

As m? is decreased further below this transition point, the stationary point of S at the 
origin i* = is seen to become progressively more unstable as the number of negative eigen- 
values of the Hessian ( p.l2| ) grows. Here the stationary points relevant to the saddle-point 



integrations ( |5.9| ) are not in the neighbourhood of the origin and are stable under changes in 
171^. The only instability that occurs with the relevant stationary points in this phase is due 
to a ^2 symmetry of the action S. Recalling its definition 

n—l n—2 

S =J2 ^i^i) - J2 ^i^i+i - ^1^ - ^n-lU (5.17) 

i=l 1=1 

we see that for vanishing x and u there is a symmetry under Zi —>■ —Zi. It follows that a 
relevant stationary point at z* has a mirror image at —z* and as x and u are tuned through 
the origin the relevant stationary point can change from one to the other. In this way non- 
analytic behaviour arises in the solution of the eigenvalue distribution po, presumably of the 
form 

Po ~ \x\ (5.18) 

which is consistent with the critical behaviour oi D = 1 matrix quantum mechanics. Of course 
this detail is difficult to verify without explicit calculation. 



In addition to this standard behaviour it was argued in |T6[ that there exist a countable 
infinity of additional points in this phase where the eigenvalue distribution po exhibits sub- 
leading non-analyticities. In fact we can reproduce these results if we assume, incorrectly, 
that the origin z = is a relevant stationary point in the evaluation of po- Considering the 
vanishing of Hessian eigenvalues (|5.12|) , we see that the naive saddle-point approximation to 
K^^^ (|5.9|) will break down when 



2 rw 



m = cos (5.19) 

s + 1 ^ ^ 
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It is well known that at degenerate stationary points there is a change in the analytic structure 
of a saddle-point approximated integral so it is natural to assume that when ( p. 19 ) is satisfied, 



K^^^ and the eigenvalue distribution po will exhibit non-analyticities. In terms of lattice 
spacing the degeneracy condition ( |5.19| ) can be found by comparing with (|5.14|) . The result 



1 - cos TTT rvr 



S + 1 



„_ „ sin— (5.20) 

V 2 2(s + l) ^ ' 

with r and s ranging over positive integers is exactly the condition for subleading critical 



behaviour as found in [|T6|. Again we stress that this behaviour is not realized in the matrix 
chain and only arises when one takes into account unstable stationary points in the evaluation 
of saddle-point integrations. This is the main shortcoming of the functional form (|3.8| ): in 
general there are a countable infinity of solutions and it is a difficult task to discover the true 
solution with minimum action. 



There is one more phase transition in the matrix chain as m^ decreases further. From (|5.12|) 



we see that as rn^ ^ — 1 all eigenvalues of the Hessian at the origin become negative and this 
stationary point is degenerate in all directions. On the surface, besides this observation of 
an irrelevant stationary point, there is little change in the system. The relevant stationary 
points are found by solving polynomial equations and their stability is determined locally as 
a smooth function of m^. This is misleading though since m^ = — 1 corresponds to a lattice 
spacing of unity which, in our scaling, is where a Kosterlitz-Thouless phase transition was 



shown to occur |]T3|, O, IIH|. The only overt sign of the transition here is that the integral 



kernel K again changes analytic structure to 

K{x,u) = e 4e 4f (5.21) 

vine 

Near the origin it has the form of a translation kernel which is of a different character than 

the original kernel ( ^.1| ). If we consider the higher, quartic terms to only make the kernel 

well-defined, the behaviour of the eigenvalue distribution pQ can be found by restricting the 

range of integration in the integral equation (|3.4| ) 



^-f /^d, p^H.) = 1 /' du e^^ e-lT'^'' ''--(^) (5.22) 

In this form it is clear that the solution will depend only on the constant b which sets a scale 



T6| , otherwise it is a universal quantity. Unfortunately the asymptotic (e -^ 0) solutions of 
this equation are not known and so remains the eigenvalue distribution pkt of the matrix 
chain at the Kosterlitz-Thouless transition. 

Finally, for m"^ < — 1 there persists a pair of !Zi2 related stationary points that at u and x 
vanishing trade off relevance. It is conjectured that this non-analytic behaviour is what one 
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would find in a D = single matrix model, i.e. 

Po ~ kP (5.23) 

It would seem that the same mechanism which produced the critical exponents of matrix 
quantum mechanics for —1 < -m^ < 1 produces different exponents here. This may be 
because of the total degeneracy of the origin but explicit calculations of the iterated kernels 



(|5.9| ) and the eigenvalue distribution po are needed to make definite statements about the 



characteristics of the solutions of this phase of the model. 

6 Conclusions and an application 

To review, we have demonstrated a general technique for expressing the equations of motion 
for systems whose time evolution is governed by Das-Jevicki-Sakita Hamiltonians with bound- 
ary conditions as an evolution problem in one-dimensional quantum mechanics. Specifying 
to the particular example of the one-dimensional chain of interacting Hermitean matrices we 
showed that in this case the general formalism reduces to a linear integral equation for the 
large- A^ eigenvalue density of the matrices. Analyzing the associated integral kernel we were 
able to recover all previous explicit solutions of the matrix chain from asymptotic solution 
of particular second order ordinary differential equations. In the general case, the solution 
of the integral equation was developed in terms of iterations of the kernel and from this 
solution the critical behavior of the system was found and a universal integral equation for 
the distribution of eigenvalues at the Kosterlitz-Thouless point was presented. In addition it 
was demonstrated that the sub-leading critical behaviour observed in [1^ results from con- 



tributions from unstable stationary points in saddle-point approximations and hence is not 
observed. 

What was not accomplished here is the explicit calculation of a non-trivial solution to the 
large- iV matrix chain using the formalism developed here. It would be nice to find a potential 
for which the series of iterated kernels could be explicitly calculated and summed and the 
eigenvalue distribution extracted. This is of particular interest in testing the conjectures 
involving the existence of a Kosterlitz-Thouless phase transition in the system. With a global 
solution for pxr one would be in a position to calculate explicitly the free energy of the matrix 
chain and check the critical behaviour. Of course only the leading non-analytic part of the free 
energy is required but in order to obtain this one should have a definite idea of the analytic 
structure of the distribution of eigenvalues. 

The calculations we have performed here are not only relevant to string theory in one 
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dimension but also for other problems related to the matrix chain. An example is the Kazakov- 
Migdal lattice model of induced QCD in D-dimensions [0]. Devised as a lattice model of gauge 
interactions where the self-interaction of adjoint scalar fields 0^ are supposed to induce the 
standard Wilson term which is left out, it is defined by the following action 

SkM = Yl ^1^(0x0 - Yj ^^^xy(py^yx(px (6.1) 

X {xy) 

Here the sub-scripts label sites on a hyper-cubic lattice on which the N x N adjoint scalar 
fields (px reside. Nearest neighbour interactions are mediated by SU{N) gauge fields A residing 
on the lattice links. From the large N saddle-point equations it can be shown |jlO| that the 
one-dimensional matrix chain we have considered ( p.l|) , with potential U, is equivalent to the 
Kazakov-Migdal model with a potential satisfying 

V'(x) = DU'ix) - 2(D - l)-fdz^^ (6.2) 

J X — z 

where the eigenvalue density po is common to both models. 

The one drawback of this model of induced QCD is that the action (|6.1| ) has a 7Zin 



symmetry under the center of the gauge group which leads to super-confinement of gauge 
degrees of freedom, even on the lattice scale |2^ . Obviously this not a feature that one would 



like to have in the continuum limit of the model and the easiest way to avoid it is for the 
system not to realize this symmetry faithfully in its solution. Unfortunately for the potentials 
U of logarithmic form, for which we have po explicitly (see Section Four), this symmetry 
persists [^, ^. It would be interesting to see if solutions of the matrix chain for more 
general potentials might lead to non-trivial propagation of gauge fields in the continuum limit 
of the Kazakov-Migdal model. This can be checked in principle using the solution of the 
matrix chain for the eigenvalue density po and the formalisms of [^ and |^ which express 
the correlation of gauge fields in the model in terms of po. 
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